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ABSTRACT 

We discuss constraints on the rates of stirring and dissipation of MHD turbulence in molecular clouds. 
Recent MHD simulations suggest that turbulence in clouds decays rapidly, thus providing a significant 
source of energy input, particularly if driven at small scales by, for example, bipolar outflows. We quantify 
the heating rates by combining the linewidth-size relations, which describe global cloud properties, 
with numerically determined dissipation rates. We argue that, if cloud turbulence is driven on small 
internal scales, the 12 CO flux (enhanced by emission from weakly supersonic shocks) will be much 
larger than observed; this, in turn, would imply excitation temperatures significantly above observed 
values. We reach two conclusions: (1) small-scale driving by bipolar outflows cannot possibly account 
for cloud support and yield long-lived clouds, unless the published MHD dissipation rates are seriously 
overestimated; (2) driving on large scales (comparable to the cloud size) is much more viable from an 
energetic standpoint, and if the actual net dissipation rate is only slightly lower than what current 
MHD simulations estimate, then the observationally inferred lifetimes and apparent virial equilibrium 
of molecular clouds can be explained. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics - ISM: magnetic fields - MHD - 
turbulence - waves 
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1. INTRODUCTION 

Millimeter-wave surveys of CO line emission have estab- 
lished the large-scale distribution and properties of molec- 
ular clouds in the Galaxy (Solomon, Sanders, & Scoville 
1979; Leung, Kutner, & Mead 1982; Dame et al. 1986; 
Solomon et al. 1987). These surveys reveal two essen- 
tial features of clouds: (1) supersonic linewidths, with 
observed velocity dispersion u exceeding the isothermal 
sound speed (c s « 0.2 km s _1 , corresponding to the mea- 
sured temperature T « 10 K) by factors ~ 10 or more; 
and (2) power-law relations between a and the cloud size 
L, of the approximate form a oc i 5 , and between the 
mean density n and L, of the approximate form n oc L^ 1 
(hereafter collectively labeled the linewidth-size- density re- 
lations). The magnitude of the observed linewidths is es- 
sentially that required to put individual clouds in virial 
equilibrium, given the estimated sizes and masses (e.g., 
Larson 1981; Myers 1983; Solomon et al. 1987). 

Turbulent motions are often invoked as an origin for 
the supersonic linewidths since they can provide an ef- 
fective pressure (Chandrasekhar 1951) which, unlike the 
much lower thermal pressure, is adequate to provide cloud 
support. However, supersonic hydrodynamic turbulence 
should decay rapidly through shocks (e.g., Mestel 1965; 
Goldreich & Kwan 1974). Since molecular clouds tend to 
have large scale magnetic fields (see Crutcher 1999), an 
attractive possibility has been that the turbulence repre- 
sents supersonic but sub-Alfvenic magnetohydrodynamic 
(MHD) waves, of which the non-compressive shear Alfven 
mode may be a long-lived component (Arons & Max 1975; 
Mouschovias 1975). 



However, recent numerical MHD simulations suggest 
that a spectrum of MHD waves, or MHD turbulence, can- 
not persist in a cloud for many crossing times (e.g., Stone, 
Gammie, & Ostriker 1998; Mac Low et al. 1998). These 
simulations tend to show that compressive motions and 
shocks arise so readily that the turbulence decays in one 
(or even much less than one) turbulent crossing time of 
the cloud (t c = L/a). This result is particularly surpris- 
ing in the simulations of Stone et al. (1998), who stir the 
system only through incompressible modes (V ■ Sv = 0). If 
confirmed, these results imply that molecular cloud turbu- 
lence must be constantly and strongly driven if clouds sur- 
vive for up to several crossing times, as implied by the in- 
ferred virial equilibrium of the clouds and their estimated 
lifetimes (e.g., Blitz & Shu 1980; Williams & McKee 1997). 

But are such high rates of dissipation really allowed, 
even if a steady-state is maintained by equally strong stir- 
ring? In this paper, we consider this question by exam- 
ining the energetics of clouds with turbulent dissipation. 
In § 2, we determine heating rates in clouds obeying the 
lincwidth-size-density relations based on dissipation rates 
determined from numerical simulations. In § 3, the corre- 
sponding luminosities are determined and compared with 
observed luminosities. Further constraints on stirring are 
given in § 4, and our conclusions are presented in § 5. 



2. HEATING RATES 

Define T(r) as the heating rate per unit volume as a 
function of position r in a molecular cloud. Various pro- 
cesses can contribute to T(r). Thermal energy input from 
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cosmic rays at the mean rate 



6.4 x 1CT 25 (n/10 3 cirr 3 ) erg cm" 3 s" 1 



(1) 



has been thought to maintain clouds of mean number den- 
sity n at temperatures T ~ 10 K (e.g., Goldsmith & Langer 
1978). 

A second contribution could also come from the decay of 
the MHD turbulence that is thought to support molecular 
clouds globally. The turbulent heating rate r tm -b can be 
written as the ratio of the energy density U of the cloud to 
an arbitrary dissipation timescale td- r tur b = U/td- The 
results of recent numerical simulations (Stone et al. 1998; 
Mac Low 1999) suggest that the ratio of dissipation time 
to turbulent crossing time k = td/t c < 1 under typical 
conditions in molecular clouds. 

For a cloud with mean density p and one-dimensional 
non-thermal velocity dispersion er, the turbulent energy 
density is U = 3/2/xr 2 , so that the turbulent heating rate 



turb 



3pa^ 
2 td 



L 



(2) 



There are a variety of other processes which can con- 
tribute to the energy balance in molecular clouds. They 
include the Galactic ultraviolet background which heats 
the edges of clouds and the internal radiation field gener- 
ated by the stars which form in localized regions of clouds. 
Here we neglect processes which might contribute only to 
localized regions of star formation and cloud boundaries. 

2.1. Stirring Scale 

An important property of the turbulent dissipation rate 
is its dependence on stirring scale. Dimcnsionally, the dis- 
sipation rate 



Tturb = V ■ 



per 
A 



(3) 



where A is the driving scale and f] is, in general, a dimen- 
sionless function of p, a and A. Kolmogorov and Obukhov 
(e.g., Landau & Lifshitz 1987) originally derived the value 
of rj in a self-similar cascade which is appropriate to the 
decay of incompressible hydrodynamic turbulence. 

By analogy, Stone et al. (1998) and Mac Low (1999) de- 
rive values of rj for their simulations of compressible MHD 
turbulence. In this case, p is the mean density of the com- 
putational region. Their results indicate that rj ~ 1 over 
a range of driving scales 1 , implying that td is of order the 
crossing time t c or less. In particular, comparing equations 
(§) and (|) shows that 
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(4) 



2.2. Dissipation in Clouds 

The linewidth-size relation is an empirically derived 
scaling between the velocity width of 12 CO ( J = 1 — 0) 
lines observed in a particular molecular cloud and the size 
of that cloud. Clouds approximately obey the scaling de- 
termined by Solomon et al. (1987): 



a = 0.72 (i?/pc) 5 km s 



-l 



(5) 



where R « is an effective radius of the cloud. The 
same study also established an approximate density-size 
relation 

p = 134 (.R/pc) -1 M pc" 3 . (6) 

These relations imply that the most massive molecular 
clouds are approximately gravitationally bound and in 
virial equilibrium. Although clouds have very complicated 
density structures, these relations give a reasonable esti- 
mate of the mean density, mean turbulent energy, and size 
of a cloud. Thus, they allow an estimate of the bolomet- 
ric luminosity of a cloud, given a mean dissipation rate as 
determined, say, from numerical simulations. 

By combining the linewidth-size and density-size rela- 
tions, we obtain the mean heating rate due to dissipation 
of energy in turbulence: 

r t urb = 7.1 x 10~ n erg cm" 3 s" 1 (7) 
= 5.4 x 10~ 25 kT 1 (ti,/10 3 cuT 3 ) 1/2 ergcm _3 s _1 . 

This relation should hold roughly for clouds with mean 
density in the range 10 2 cm~ 3 < n < 10 4 cm -3 . Also note 
that the energy density (or pressure) is independent of 
cloud size. 

3. CLOUD LUMINOSITIES 

Assume only that the cloud exists in some definite vol- 
ume and that it is in a steady-state equilibrium. Under 
these circumstances, the total cooling rate within the vol- 
ume must balance the total heating rate. Therefore, al- 
though the cloud may have a rather complex density struc- 
ture, and a variety of atomic and molecular species (each 
having a different optical depth) may contribute to the 
cooling, the bolometric luminosity Lboi must balance the 
volume-integrated heating rate: 



U = JdrT(r)=L boh 



(8) 



Note that we have neglected the possibility that mass loss 
can carry away excess kinetic energy from the cloud in the 
form of a global wind or outflow. 

Of course, determining the bolometric luminosity of a 
cloud is no trivial task. Current observations usually only 
consider specific atomic and molecular line transitions, 
with 12 CO (J = 1 — 0) by far the most thoroughly in- 
vestigated. In their study, Solomon et al. (1987) deter- 
mine a 12 CO (J = 1 — 0) luminosity- linewidth relation 



Since driving may occur on scales A much smaller than 
typical cloud sizes L ~ few — 100 pc (e.g., through out- 
flows), the results of numerical simulations suggest that k 
may be much less than unity. This can obviously result 
in very large cloud luminosities; below we determine the 
expected fluxes for these heating rates. 

1 Mac Low (1999) defines a constant r\ v using the three-dimensional velocity dispersion and other numerical factors, so that our rj = 32.6r/ v 



for the clouds in their sample (their eq. [9]). We use the 
linewidth-size-density relations and convert units in their 
formula to express the empirical 12 CO ( J = 1 — 0) lumi- 
nosity as 



CO 



3.1 x 10 30 (n/10 3 cm~ 3 )- 5/2 



ergs 



(9) 
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This luminosity represents only the flux through the por- 
tion of the cloud's surface which is visible to the observer. 
Here we adopt a total CO luminosity Leo = 4L' co , ap- 
propriate for a spherical cloud. Of course, real molecular 
clouds have irregular shapes and surfaces, so that the mul- 
tiplicative constant may be somewhat greater than 4, but 
it should be no more than an order unity correction. 

For clouds whose average properties obey the linewidth- 
size-density relations, the bolometric luminosity deter- 
mined from equation d8h is 



Lboi = 8.5x10 



:',2 



1.2 (n/10 3 )- 2 + kT 1 (n/10 3 )~ 5/2 



ergs 
(10) 

where the first term in brackets gives the cosmic ray heat- 
ing from equation (|l|) and the second term gives the tur- 
bulent heating from equation (Q). 
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Fig. 1 Top: Bolometric luminosity Lboi (equal to the sum 
of cosmic ray heating and turbulent dissipation) as a func- 
tion of mean density n for clouds which obey the linewidth- 
size-density relation and have fixed values of the dissipa- 
tion parameter re (dashed lines) or the turbulent driving 
scale A (dotted lines). The solid line shows the 12 CO lu- 
minosity Leo from the empirical relation determined by 
Solomon et al. (1987). Bottom: Ratio of observed Leo 
to bolometric luminosity Lboi for the different values of re 
and A. 

The top panel of figure 1 shows the bolometric lumi- 
nosity (equivalent to the heating rate) compared with the 
expected Leo from Solomon et al. (1987). Dashed lines 
show Lboi for two values of re, and dotted lines show Lboi 
for two values of a fixed driving scale A that are compatible 
with driving by bipolar outflows. The bottom panel shows 
the flux ratio Leo /Lboi, thereby illustrating how much 



emission is required beyond the 12 CO J — 1 — transition 
to balance the heating rate. The excess can reach (or even 
exceed) 3 — 4 orders of magnitude, depending on the cloud 
mean density and driving scale. Detailed models of molec- 
ular line cooling (Neufeld, Lepp, & Melnick 1995) suggest 
that Lco/Lboi typically falls in the range 1% — 2%. There- 
fore, it may be difficult to account for an additional 1 — 2 
orders of magnitude excess emission. Yet this is required 
if clouds are stirred from within, since bipolar outflows 
are usually invoked as the only viable stirring mechanism 
(Norman & Silk 1980; McKee 1989; Shu et al. 1999), and 
these act on relatively small scales A w 0.1 — 1 pc. 

Can the required excess emission be hidden in the high 
J transitions of CO, or in lines from other species? In the 
recent simulations of Stone et al. (1998), Smith, Mac Low, 
& Zuev (2000; hereafter SMZ), and Smith, Mac Low, & 
Heitsch (2000; hereafter SMH), some 50% -68% of the tur- 
bulent energy dissipates in shocks. In simulations of freely 
decaying MHD turbulence, SMZ find that most shocks are 
slow, with typical Mach number in the range 1 — 3, i.e., 
velocities in the range v s = 0.2 — 0.6 km s _1 for a kinetic 
temperature T — 10 K. Applying the Rankine-Hugoniot 
jump conditions for a nonradiative shock, this range of v s 
yields post-shock temperatures at best as high as T — 27 K 
(for v s — 0.6 km s _1 ). SMH find only slightly higher 
Mach number velocity jumps for driven MHD turbulence 
in systems with an energetically significant magnetic field. 
Even if shock velocities are as high as v s — 1 km s _1 , 
the maximum possible (nonradiative) post-shock temper- 
ature is T = 87 K, which lies within the range in which 
emission from CO dominates the molecular emission for 
n < 10 4 cm -3 (see Fig. 4 of Neufeld et al. 1995). Further- 
more, the low order transitions of CO are easily excited 
under these conditions, with the lowest transition of 12 CO 
(J = 1 — 0) accounting for Leo ~ 10~ 2 Lb o i- (Neufeld 
et al. 1995). Therefore, the high bolometric luminosities 
Lboi shown in Figure 1 for small scale driving (implying an 
observed Lco/Lboi ~ 10~ 4 — 10~ 3 ) are incompatible with 
molecular cooling models. It is likely that a detailed cal- 
culation based on the distribution of shock strengths in a 
typical simulation with small-scale driving will show a sig- 
nificant excess of Leo above the empirical values presented 
in Figure 1. This will also imply much greater values of 
the excitation temperature than actually observed. 

Finally, we note the scalings of the luminosity and en- 
ergy generation rate. Since cloud luminosities scale as 
n~ 2 - 5 (alternatively L 2 5 , leading to the mass-luminosity 
scaling M oc L 8 measured by Solomon et al. 1987), 
it may be difficult to attribute their radiative output to 
cosmic ray heating, which on average scales as n~ 2 (or 
L 2 ). Turbulent dissipation proportional to the crossing 
time provides the correct scaling for energy input, i.e., 
when the driving scale is proportional to the individual 
cloud size L rather than a fixed internal scale A < L. In 
fact, dissipation with re J> 1, consistent with A sa L, yields 

Lco/Lboi ~ 10 -2 (see Fig. 1), roughly compatible with 
the observed Leo and the application of molecular cool- 
ing models to post-shock conditions of weakly supersonic 
shocks. 



4. OTHER CONSTRAINTS 
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In addition to the luminosity problem, another feature 
of rapid internal driving is that an amount of energy in ex- 
cess of the gravitational binding energy of a cloud needs to 
be input and dissipated every crossing time (since A <C L). 
This is not likely to be a very stable situation, as any frac- 
tional imbalance between dissipation and cooling can lead 
to expansion and mass flow from the cloud. 

Furthermore, the star formation rate required to sup- 
port this energy input also appears to be too large. An 
estimate for the required star formation rate can be ob- 
tained in the manner first done by McKee (1989) (see also 
Shu et al. 1999), but incorporating the published dissi- 
pation rates of Stone et al. (1998) and MacLow (1999), 
and using the linewidth-size-density relations. The local 
dissipation rate given by equation (||), with 77 = 1, yields 
a Galactic dissipation rate 



L 



turb 



Ma 3 
A 



(11) 



where M is the mass of molecular gas in the Galaxy. We 
compare this to the input rate from momentum-driven 
bipolar outflows with launching speed v w and mass flux 
M w which equals a fraction / of the star formation accre- 
tion rate Msf- If the swept-up shell effectively dissipates 
when the speed drops to i; rms = V3a, the delivered energy 
is 

Lout = -^fM SF v w v rms , (12) 

Equating equations ( |TT| ) and (|l2|), the required star for- 
mation rate is 



M SF = 



Ma 2 



3V2/ \ Vv 



(13) 



To obtain a numerical estimate, we use / = 1/2 (e.g., 
Shu et al. 1999), v w — 200 km s _1 , and a relatively 
large driving scale A = 1 pc for the outflows. The to- 
tal mass M ~ 10 9 Af Q of molecular gas in the Galaxy is 
concentrated in the giant molecular clouds (e.g., Solomon 
et al. 1987), for which we take n = 10 2 cm~ 3 to be a 
representative density. Equations (||) and @ then yield 
a = 3.48 km s -1 . Finally, evaluating equation (|l^) yields 
Msf — 143M Q yr _1 , which is clearly too high compared 
to most estimates of the Galactic star formation rate 
3- SMoyr" 1 (see McKee 1989). Altogether, recov- 
ering a low enough star formation rate given the fast dis- 
sipation rates can be done only with extreme fine tuning, 
and therefore seems unlikely. Our conclusion is similar to 
that of McKee (1989) that bipolar driving is adequate if 
the dissipation occurs over many free-fall times tg of the 
cloud, but is inadequate if it occurs over ~ ltg . 

5. CONCLUSIONS 

In molecular clouds, energy generation by turbulent 
dissipation can far exceed heating by cosmic rays (c.f. 
eq. ]Ic|]), which have traditionally been viewed as the 



energy source responsible for the 10 K excitation tem- 
peratures measured in 12 CO (Goldsmith & Langer 1978; 
Neufeld et al. 1995). In this paper, we have demonstrated 
that rapid internal turbulent driving on small scales, in 
clouds which obey the linewidth-size-density relation, and 
in which MHD turbulence damps on essentially a cross- 
ing time over the driving scale, requires that the observed 
12 CO ( J = 1 — 0) luminosity Lqq undersample the bolo- 
metric cloud luminosity Lboi by a far greater factor than 
predicted in standard molecular cooling models (Neufeld 
et al. 1995). We are left with two distinct possibilities: 
either current numerical simulations are seriously overes- 
timating the dissipation rate of MHD turbulence, or the 
idea of long-lived molecular clouds supported for several 
global crossing times by small-scale internal driving, is in- 
correct. 

If stirring occurs on approximately the scale L of each 
individual cloud, then the observed magnitude and scal- 
ing of Leo can be understood with MHD dissipation rates 
comparable to those calculated in recent simulations. This 
may mean that clouds are not long-lived and need not 
maintain a steady-state through stirring equal to the dis- 
sipation rate, consistent with the conclusion (based on stel- 
lar ages) that star formation occurs on a single cloud cross- 
ing time (Elmegreen 2000). Alternatively, clouds need to 
survive for only « 2 — 3 crossing times to explain their 
estimated lifetimes and their apparent state of virial equi- 
librium. We note that the crossing time for clouds which 
obey the linewidth-size-density relation is 



t c = L/a = 1.3 x 10 7 (n/10 2 cm- 3 )- 1/2 



which is ~ 4 times the free-fall time 



vr, 



(14) 



< ff = (37r/32G/9) 1/2 = 3.3 x 10 6 (n/10 2 cm- 3 )-^ 2 yr. (15) 

The calculated t c is only a factor « 2 — 3 lower than the 
estimated lifetime w 3 x 10 7 yr of giant molecular clouds 
(Williams & McKee 1997) of mean density n w 10 2 cm~ 3 . 
Thus, we believe that future numerical simulations which 
yield slightly (a factor w 2 — 3) lower net dissipation rates 
can reconcile turbulent dissipation (if driving occurs on 
the scale L of each cloud) with the apparent evidence for 
long-lived (lifetime « 2 — 3< c ~ 10 tg) clouds that are in 
near virial equilibrium. These somewhat lower dissipation 
rates may be obtained by tapping previously unmodeled 
global sources of turbulence such as gravitational contrac- 
tion and participation in the differential rotation of the 
Galaxy, and/or dealing more carefully with heating and 
cooling of low density molecular cloud envelopes. The lat- 
ter requires dropping the isothermal assumption, which 
can overestimate radiative losses and enhance the devel- 
opment of shocks, particularly in low density regions. 

We thank Ted Bergin for discussions and encourage- 
ment, and an anonymous referee for insightful comments. 
SB was supported by a grant from the Natural Sciences 
and Engineering Research Council of Canada. 
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